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Abstract. A single fermion is bound by a softcore central Coulomb potential 

V{r) = — D/(r' +b'^)i , V > 0, b > 0, q > 1, in d > 1 spatial dimensions. 
Envelope theory is used to construct analytic lower bounds for the discrete Dirac 
energy spectrum. The results are compared to accurate eigenvalues obtained 
numerically. 
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1. Introduction 

For a central potential in d dimensions the Dirac equation can be written [T in natural 
units h = c — 1 as 

i— — = H'^, where H = ctsPs + nT-P + V, 

OT 

s — 1 

m is the mass of the particle, ^ is a spherically symmetric vector potential, and 
{oLs\ and /? are the Dirac matrices which satisfy anti-commutation relations; the 
identity matrix is implied after the potential V . For stationary states, some algebraic 
calculations in a suitable basis lead to a pair of first-order linear differential equations 
in two radial functions {"01 ("r), V'2 ('')}, where r = ||r||. For d > 1, these functions 
vanish at r = 0, and, for bound states, they may be normalized by the relation 

oo 

(V'l, Vi) + (^2>2) = J ii^lir) + ^l{r))dr = 1. 



We use inner products without the radial measure ' because the factor r 2 is 
already built into each radial function. It has been shown, for example by Jiang [T] 
(using algebraic ladder-operator methods), that these functions satisfy the following 
coupled radial equations 

E^Pi = {V + m)iJi + {-d + kd/r)4,2, (1) 
E^2 - (5 + kd/r)^Ji + {V-m)^P2, (2) 

where d = d/dr, ki = 0, fc^ = T{j + d > 1, r ±1, and j = 1/2,3/2,5/2,... 

is the total angular momentum quantum number. We note that the variable r is 
sometimes written as w, as, for example in the book by Messiah [2], and the radial 
functions are often written as ^"1 = G and '02 = F, as in the book by Greiner [3]. 
We shall assume that the potential V is such that there is a discrete eigenvalue 
E and that Eqs.([l] [2]) are the eigen-equations for the corresponding radial eigen- 
states. Our geometrical method (to be outlined below) will also presume that the 
potential can be written as a smooth transformation V{r) — g{—l/r), where g is 
monotonically increasing and of definite convexity. In this paper we shall present the 
problem explicitly for the cases d > I. We shall label the discrete eigenvalues by kd 
and the number = i^i of nodes in the large radial component ipi. This convenient 
labelling is suggested by a result of Rose and Newton [4j to the effect that if 1/2 is the 
number of nodes in the small component "02, then if r = —1, i>2 = t^i, and if r = 1, 
then 1^2 — 1^1 + ^- We study the discrete Dirac spectrum generated by the potential 
V{r) = Vq{r) = vfq{r), where ?; > is the coupling parameter, and the potential 
shape fq(r) is given by 

This potential represents a family of softcore (truncated) Coulomb potentials, which 
are useful as model potentials in atomic and molecular physics. In the limit as q ~> 00, 
the potential descends to the cut-off Coulomb potential f^o given by 



lim fg{r) = foo{r) 

q^oo 



i, ifr<6; 



ifr > b. 
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The bound states are obtained in terms of three potential parameters: the 
coupUng > 0, the cut-off parameter 6 > 0, and the power parameter q > I. The cases 
g = 1 and 2 are of special physical significance [51-[T5]. The potential /i represents 
the potential due to a smeared charge and is useful in describing mesonic atoms. 
The potential /2 is similar to the shape of the potential due to a finite nucleus and 
experienced by the muon in a muonic atom. Extensive applications of the softcore 
Coulomb potential, /2, have been made through model calculations corresponding to 
the interaction of intense laser fields with atoms p^HlQ] . The parameter b can be 
related to the strength of the laser field, with the range 6 = 20 — 40 covering the 
experimental laser field strengths [14 . 

In the non-relativistic case, Mehta and Patil [5] have presented analytical solutions 
for the s-state eigenvalues corresponding to the /i potential. Also upper and lower 
energy bounds and some exact analytic solutions have been found for special cases 
with the potentials /i and /2 [201 121] • Patil [B] has discussed the analyticity of the 
scattering phase shifts for two particles interacting through the potentials fq with q = 1 
and q = 2. Much less is known concerning the corresponding relativistic problem. 

The principal idea that is used in this paper is that of envelope theory. We suppose 
that a given potential shape /(r) can be represented as a smooth transformation 
/(r) = g{h{r)), where h{r) is a potential that generates a soluble spectral problem. 
Since tangents to g{h) are of the form ah{r) + c, they generate a family of soluble 
'tangential problems'. If the transformation function g has definite convexity, these 
tangential problems lead via comparison theorems to a set of energy bounds. Envelope 
theory picks out the best of these. Although the principal focus of the present paper 
is on spatial dimensions d > 1, we note that in d 1 dimension, the corresponding 
potential would be of the form V{x) = vf{\x\), and the spatial components ^i{x) 
and ip2{x) of the spinor can be classified as even or odd functions. These components 
need not now have to vanish at the origin unless the potential is sufficiently singular 
there. Meanwhile, the normalization would be given by J_^{ipiix) + V'K^^))'^^; = 1. 
The geometrical reasoning used in the present paper applies equally well to problems 
in one dimension but, of course, necessitates useful exactly soluble problems to be 
used as bases for the approximations; the Dirac Coulomb problem, for example, is 
problematic in d = 1 dimension |22 l f23 ] . Nieto [2^ has presented an analysis of non- 
relativistic problems in d dimensions, with d real and positive. For the future, the 
task of discussing a similar generalization presents itself for relativistic problems, and 
smooth transformations thereof. 

In section 2 we review some recent comparison theorems for Dirac eigenvalues, 
and we discuss general scaling and monotonicity properties. In section 3 we describe 
envelope theory for the Dirac equation, and in section 4 we look at the specific case in 
which the softcore Coulomb potential is written as a smooth convex transformation of 
a pure Coulomb potential: this generates a simple formula for lower energy bounds. 
We find some of these bounds explicitly and compare them with accurate values found 
by direct numerical methods. 

2. Comparisons, monotonicity, and scaling 

It is not a simple matter to characterize the discrete Dirac spectrum variationally 
|25fl27) . However, in spite of this, some comparison theorems have recently been 
proved [28H30] , and we state two of these theorems here for use in the present paper. 
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Theorem 1 |28| The real attractive central potential V{r, a) depends smoothly on the 
parameter a, and E{a) — Ek^j^{a) is a corresponding discrete Dirac eigenvalue. Then: 

dV/da > ^ E'{a) > and dV/da < ^ E'{a) < 0. 

Theorem 2 |29] Suppose that E^^^^ and E^J^ are Dirac eigenvalues corresponding to 
two distinct attractive central potentials V''^^{r) and V^^\r). Then: 

V^'\r)<V^^\r) <i<<l. 

If the exact eigenvalues of H are written E(v,b,q,m), then we conclude from 
Theorem 1 and the monotone behavior of the potential V{r) = vfq{r) with respect 
to the parameters that these spectral functions are monotone in each parameter, 
decreasing in v and q, and increasing in b. Because 

dE ^ dE ^ , dE ^ 
^ < 0' 1^ > 0, and ^ < 0. 
ov ob oq 

We now change variable r ^ Srva Eqs.(ITJ[5]), where (5 > is constant, multiply through 
by 5, and compare eigenvalues, we obtain the general scaling law for the family of 
softcore Coulomb potentials V{r) = vfq(r) under the Dirac coupled equations, namely 

E{v, b, q, m) = -E ^u, -, q, 5m 

Choosing 6 — b and 6 ^ 1/m we get, respectively, the special scaling laws 

E(v, b, q, m) — —E (v, 1, q, bm) = mE (v, bm, q, 1) . 



This much is known in general for the whole class of problems. 



3. Envelope theory 

Envelope theory has been used since 1980 as a geometrical method of spectral 
approximation [31-35 . We include here a brief self-contained summary of what is 
needed for the present task. We consider a potential shape /(r) that can be written 
as a smooth transformation /(r) — g{h{r)) of a potential h{r) for which the solutions of 
the Dirac equation are exactly known. In our specific application, h{r) will be chosen 
as the Coulomb potential h{r) = — l/r, and /(r) = ,fq{r) is the softcore Coulomb 
potential, which we are studying. However, it is clearer to discuss the method in 
general at first. Let us suppose that the transformation function g{h) is monotonically 
increasing and convex, that is to say, g'{h) > and g"{h) > 0. This means that g{h) 
lies above its tangents. If r = i is the point of contact between curve and tangent, we 
have (in the convex case) a family of 'lower' tangential potentials given by 

f{r)>f('\r)^a{t)h{r)+c{t), (3) 

where 

ait) = and c{t) ^ g{h{t)) - h{t)g' {h{t)). (4) 

We note parenthetically that if g{h) is concave^ we obtain, instead, a family of 
upper bounds. Continuing the convex case, \i E = D(u) describes how a discrete 
Dirac eigenvalue corresponding to the potential uhir) depends on the coupling u, 
then the potential inequality Eq. Q and Theorem 2 imply for the original potential 
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V{r) — vf{r) the spectral inequality E > D{va{t)) + vc{t). This expression can then 
be optimized over the contact point t to give the lower bound 

E{v) > E^{v) = sup [D{va{t)) + vc{t)] . (5) 

Replacement a{t) and c(t), using Eq. in Eq. ([5]) and differentiation with respect to 
h gives the value for the critical point h = D'{vg'{h)). Then we can re- write the right 
hand side of Eq. ([5]) by changing the minimization variable from t to u by means of 
the invertible transformation u = vg'{h(t)), yielding the following alternative form for 
the best energy bound: 

E{v) > E^{v) = sup [D{u) - uD'{u) + vg{D'{u))] . (6) 
«>o 



4. Energy bounds for the softcore Coulomb potential 



We first consider the Dirac equation for the pure Coulomb problem with potential 
V{r) = —u/r^ where the coupling parameter u = aZ is not too large. We write the 
exact discrete eigenvalues as Dk^v{u) = D{u) and they are given pQ[22] exactly by 



(7) 



D{u) = m\l + u-' v-{l-T)/2 + {k%- u^)" 



where 



d-2 



n-|fcd| + (fc^-u^)3 



< u < 1, 



t[j + 



±1, 



= 0, 1, 2, ... is the number of nodes in the upper radial function 'ipi{r), and n is the 
principal quantum number defined in general (for both Coulomb and non-Coulomb 
central potentials) by 

n = iy+ \kd\ —■ 

The spectroscopic designation 

{s,p,d,...} ^£^{0,1,2,...} 
is then provided by the formula 

For the softcore Coulomb potential fq{r) it follows that the transformation 
function g{h) = fq{—\/h) has the explicit form 

9{h) = -— ^ -T. (8) 

Consequently dg/dh and d^g/dh^ are given by 



dg 



dh h^i-iY + b^] 



1/9+1 



and 



d'g 



1/9+2 ■ 
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Figure 1. The softcore Coulomb potential for q = 2 shown with tangential family 
of shifted Coulomb potentials. 

Since the quantity —1/h > 0, we conclude that dg/dh > and d^g/dh^ > 0. Thus g 
is monotonically increasing and convex; its tangents f^*\r) = a{t)h{r) + h{t) given by 
Eq. dH) are lower bounds to fq{r). The case g = 2, = 2, and w = 58 is shown in Fig. 
©■ 

Lower energy bounds are therefore provided by Eq. ([6|) with D(u) given by Eq. ([7]) 
and the transformation function g given explicitly by Eq. ([8|) . The lower bounds ^ 
are compared with accurate numerical values Ek^u in Table 1. 

The eigenvalue formula Eq. ([7]) is based on the Coulomb spectral function D{u) 
which admits couplings u satisfying u < 1. Thus we require in Eq. ([S]), that 
va{t) < 1 and consequently we cannot consider arbitrarily large coupling v for the 
softcore Coulomb potential. The Coulomb degeneracy -Efe^,/ = ^'-fc^iy+i is expressed 
by Eq. ([7]) and, of course, this symmetry is satisfied by the lower bounds. For 
instance, from Table 1 we see that the lower bounds E^^^ are degenerate for the 
pairs {i?i33, -E_i34} and {-E53I, £'-532}, although the eigenvalues themselves are not 
exactly equal. Nevertheless, the simple lower bound formula Eq. (|6]) is valid for all the 
discrete Dirac eigenvalues and is often very effective. 
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Table 1. Comparison the exact softcore Coulomb energy eigenvalues E with the 
lower bounds , for m = 1. 
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5. Conclusion 

This paper is based on two ideas: (1) a comparison theorem vaUd for a discrete 

Dirac spectrum, and (2) a geometrical theory that generates a family of potentials 
tangential to a smooth transformation g{h) of a Dirac-soluble base potential h. These 
two strands are connected if g{h) has definite convexity so that its graph lies either 
above or below its tangents. Meanwhile the tangential potentials are of the form 
ah{r) + c and Dirac's equation is soluble exactly for each of them. For given values 
of the potential parameters, and good quantum numbers, one is then able to find the 
best tangential potential in the sense of providing the best energy bound. We have 
used the Coulomb envelope base h(r) = — 1/r to generate optimized energy lower 
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bounds for a family of softcore Coulomb potentials given by g{h{r)) = —l/{r'' + 
where 6 > and q > 1. Bccaiiso the potential is central, the geometric argument 
leading to the lower bound via the Dirac comparison theorem transcends the question 
of the number d of spatial dimensions. The estimate is an energy bound whenever g 
has definite convexity. 
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